Determining the surface-to-bulk progression in the normal-state electronic structure 
of Sr 2 Ru04 by angle-resolved photoemission and density functional theory 

C.N. Veenstra, 1 Z.-H. Zhu, 1 B. Ludbrook, 1 M. Capsoni, 1 G. Levy, 1 ' 2 A. Nicolaou, 1 ' 2 

J. A. Rosen, 1 R. Comin, 1 S. Kittaka, 3 Y. Maeno, 3 I. S. Elfimov, 1 ' 2 and A. Damascelli 1 ' 2 ' * 

1 Department of Physics & Astronomy, University of British Columbia, Vancouver, British Columbia V6T1Z1, Canada 
2 Quantum Matter Institute, University of British Columbia, Vancouver, British Columbia V6T1Z4, Canada 
3 Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan 

(Dated: February 14, 2013) 

We revisit the normal-state electronic structure of Sr2Ru04 by angle-resolved photoemission spec- 
troscopy (ARPES) with improved data quality, as well as ab-initio band structure calculations in the 
local-density approximation (LDA) with the inclusion of spin-orbit coupling (SO). We find that the 
current model of a single surface layer (y/2 x v / 2)R45° reconstruction does not explain all detected 
features. The observed depth-dependent signal degradation, together with the close quantitative 
agreement with LDA+SO slab calculations based on the surface crystal structure as determined 
by low-energy electron diffraction (LEED), reveal that - at a minimum - the subsurface layer also 
undergoes a similar although weaker reconstruction. This model accounts for all features - a key 
step in understanding the electronic structure - and indicates a surface-to-bulk progression of the 
electronic states driven by structural instabilities, with no evidence for other phases stemming from 
either topological bulk properties or the interplay between SO and the broken symmetry of the 
surface. 

PACS numbers: 74.25.Jb, 74.70.Ad, 79.60.Jv, 79.60.Bm 



Since its discovery [1], understanding the origin of su- 
perconductivity in Sr2Ru04 has received enormous at- 
tention [2] . In search of the prodromes of the elusive pair- 
ing mechanism, its normal-state electronic structure has 
also been a matter of intense study. A unified picture of 
the Fermi surface (FS) was obtained from bulk-sensitive 
de Haas- van Alphen (dHvA) measurements [3, 4] and sur- 
face sensitive angle-resolved photoemission spectroscopy 
(ARPES) [5]. This, together with local-density approx- 
imation (LDA) band structure calculations [6, 7], es- 
tablished the importance of spin-orbit coupling (SO) in 
defining the low-energy electronic structure. The inclu- 
sion of SO may also lie at the root of the bigger question 
regarding the superconducting mechanism, with sugges- 
tions that it may facilitate mixed-parity pairing [6, 8, 9]. 

Similar to the case of topological insulators the pres- 
ence of SO also raises the possibility that novel topo- 
logical states, Dirac or Rashba-type, might exist on the 
(001) cleaving-surface of Sr2Ru04 in the normal- state (in 
addition to the chiral superconductivity topological edge- 
states detected below T c on the side surfaces [10]). This 
proposal is particularly intriguing given that surface- 
specific electronic states have been observed in this mate- 
rial [5, 11]. While these have originally been interpreted 
in terms of a purely structural reconstruction of the top 
surface layer [5, 11] as resolved via low-energy electron 
diffraction (LEED) [12], a recent ARPES study detected 
an additional dispersive band [13], demanding this simple 
model be revisited in terms of a possible manifestation 
of topological bulk properties or Rashba interaction. 

Here we study the normal-state electronic structure 
of Sr2Ru04 by ARPES and density-functional theory 
(DFT). With improved data quality a multitude of new 



dispersive features can be detected, stemming from the 
folding of all FSs: reconstructed-surface's and bulk-like's 
as well. While this might be suggestive of surface ferro- 
magnetism or Rashba-type splitting, we show - through 




FIG. 1: (color online). The LDA+SO, k z = FS of Sr 2 Ru0 4 
for: (a) undistorted bulk; (b) bulk with Ru06 octahedra ro- 
tated 5.5° causing a (y/2 xy / 2)R45° reconstruction (new BZ - 
dashed; unfolded FS - solid; folded replica FS - dotted), (c) 
ARPES FS (7meV integration) from a fresh low-temperature 
cleave; (d) phenomenological FS based on (c), Fig. 2, and 3. 
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FIG. 2: (color online), (a) Sr^RuCU FS maps (grey scale), marking the location of MDC cuts at the Fermi energy showing 
clear doubling in the (b) f3 sheet, (c) its folded replica f3\ (d) a, and (e) folded a . (b-e) Data are shown in color, with red and 
blue referring to surface and bulk- like features, while fit results (see text) are shown in black (dashed lines for folded bands). 



controlled aging and LDA+SO slab calculations for the 
LEED-determined crystal structure [14] - that this dou- 
bling arises from surface and subsurface instabilities. 
This discovery bears close similarity to the recent findings 
of surface-enhanced charge-density- wave [15] and stripe 
[16] order in the underdoped high-T c cuprates, establish- 
ing structural/electronic instabilities in the near surface 
region a widespread phenomenon of complex oxides. 

ARPES experiments were performed at UBC with a 
SPECS Phoibos 150 analyzer and 21.2 eV photons from 
a monochromatized UVS300 lamp. FWHM energy and 
momentum resolutions were measured to be 17meV and 
0.01 ^ angular equivalent in Fig. 1 and 2, and 9meV and 
0.04 -7=- in Fig. 3. Samples grown by the floating-zone 
method [17] were oriented by Laue diffraction, cleaved 
in-situ, then held at temperatures between 5-10 K and 
pressures varied in the 5 — 9x 10 -11 mbar range. To ver- 
ify that the observed ARPES features are not associated 
with Ru-metal inclusions, the so-called '3K-phase' [17], 
samples with and without inclusions were studied: no 
3K-phase dependence was observed. 

Let us begin with a brief review of the previously ac- 
cepted electronic structure from bulk and reconstructed 
surface [5, 6, 11]. Fig. 1(a) shows the LDA+SO bulk FS 
derived from the Ru-^ orbitals occupied by 4 electrons, 
which is in good agreement with results from dHvA [3, 4], 
as well as ARPES on an optimally degraded surface [5]. 
As a consequence of the (a/2x\/2)R45 surface reconstruc- 
tion observed by LEED [12] and associated with the rota- 
tion of the Ru06 octahedra [Fig. 4(c)], the Brillouin Zone 
(BZ) is reduced leading to band folding; this is shown 
in Fig. 1(b) for a 5.5°-rotation LDA+SO bulk calcula- 
tion (note that primes and dotted lines will be used for 
folded FSs, and s/red for surface and 6/blue for bulk- 
like states). In Fig. 1(c) we present the ARPES FS for a 
pristine cleave of S^RuC^, measured in p-polarization at 
5 K and 6xl0 -11 mbar in less than 3 hours after cleaving, 
during which time no aging was observed. The data were 
acquired over the entire upper right quadrant without 



symmetrization, then folded in /c-space to reproduce the 
BZ as shown (the analyzer entrance slit was oriented hor- 
izontally with respect to the final image, giving greater 
angular resolution in that direction). One can recognize 
doubling and folding of many features, seemingly consis- 
tent with the overlap of bulk and reconstructed FSs from 
Fig. 1(a) and (b), as originally proposed [5, 11]. For in- 
stance, both fa and f3 s are seen, as well as both 75 and j s 
with their expected electron and hole-like topology [with 
7 S showing clearly near (0,7r) and 75 near (7r,0), due to 
varying matrix element dependence]. However, contrary 
to this explanation, the folded /?' also is clearly doubled 
giving rise to j3' h and j3' s . We will show in Fig. 2 and 3 that 
in fact all bands are doubled - both unfolded and folded 
- resulting in the phenomenological FS of Fig. 1(d). This 
cannot be explained by a single reconstructed surface- 
layer model, and may instead be consistent with the in- 
teresting proposal of SO-driven surface effects [13]. 

In Fig. 2 we confirm the doubling of features with mo- 
mentum distribution curves (MDCs) obtained by inte- 
grating a 10 meV region just above Ej?: panels (b-e) show 
the doubling of /3, j3 r , a, and a' , respectively, observed 
at the /c-locations indicated in the alignment FS of (a). 
Note that the splitting between a' s and a' h is more clearly 
resolved in the second BZ [Fig. 2(e)], perhaps due to ma- 
trix element effects and/or 2 nd -BZ higher angular resolu- 
tion. By performing a least-squares fit to two Lorentzian 
peaks with a second-order polynomial background (black 
curves), we can estimate the observed splittings as: 0.057, 
0.065, 0.027, and 0.029±0.003 f for /3, /3', a, and a'. Al- 
though it is difficult to compare these values directly, due 
to the slight difference in angles and /c-space locations be- 
tween them, we note that the splittings between P s /Pb 
and f3' s / f3' h are very close and approximately double those 
between a s /ab and a' s /a' h , which are almost identical. 

Thus far, with a doubling observed for all a and /3 
bands and in light of the bulk magnetism of many Ru- 
oxides, one might be drawn to the possibility of Rashba- 
type effects [13] or even surface ferromagnetism [12]. If 
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FIG. 3: (color online), (a) Time evolution of a 45° MDC 
cut at Ep\ its location is marked on the alignment FS in (b). 
(c,e,g) Plot of the corresponding MDC at Ef, band map of 
features on a large energy range, and zoom in energy on the 
Fermi crossings for the first 100 minutes after cleave; (d,f,h) 
the same for the final 100 minutes of acquisition. In (g,h), 
MDC maxima (black circles) and MDC-dispersion fits (blue- 
bulk/red-surface lines) are also shown. 



any of these were the case one would expect that - under 
surface degradation - the surface states' splitting would 
shrink and/or the intensity decay uniformly over time for 
all of them, eventually revealing the underlying bulk elec- 
tronic structure. In Fig. 3 we show the results of such a 
time-dependent experiment for the band dispersion along 
the (0, 0)— (7r, 7r) direction in s-polarization, where all fea- 
tures and their folded replicas are visible simultaneously, 
for a sample maintained at 9.5 K while raising the base 
pressure to about 8.5xl0 -11 mbar. 

The time-evolution of MDCs integrated over 100 min- 
utes and 4meV just above Ep is presented in Fig. 3(a): 
while some features remain at a similar level of intensity 
over the 14-hour time-span, others age and have their in- 
tensity suppressed, without any overall change in slope or 
shape. Particularly comparing Fig. 3(c,e,g) and (d,f,g), 
which present only first and final MDC and dispersion 
results, we can infer that: (i) all bands and their folded 



replicas are doubled for a total of 12 bands (as labelled), 
including the 7 bands not previously demonstrated in 
Fig. 2; (ii) based on the observed difference in degrada- 
tion rate, the electronic states do not all originate from 
the very same surface layer, and can instead be classified 
as surface (fast degradation) or bulk-like (slow degrada- 
tion) . These observations can be made more quantitative 
from the MDC-dispersion fit analysis, which allows dis- 
tinguishing between surface and bulk-like bands based on 
their remarkably different Fermi velocity vp (see Tab. I). 
While the fresh data in Fig. 3(g) are dominated by the 
surface a s /a' s , P s /P' s , and 7 s /7g, with the addition of the 
well-separated Pb/P' h pair, the aged data of Fig. 3(h) only 
show the much steeper bulk-like at,/a f b , Pb/ P' h , and 75/75 
(a r h , although visible, could not be fit and is shown with 
a dashed line as a guide to the eye). 

Since all s-states are suppressed evenly and indepen- 
dently from the 6-ones, we must conclude that neither 
ferromagnetism or Rashba coupling is responsible for the 
apparent band doubling. Similarly, we must rule out pos- 
sible surface patches with different structure or effective 
doping, as well as bulk photoelectrons being scattered off 
the reconstructed surface [18]. As intrinsic bulk proper- 
ties can also be excluded, otherwise they would be visible 
in dHvA [3, 4] and x-ray diffraction [19], the most natu- 
ral explanation would be that the subsurface layer has a 
structural modulation similar to, yet weaker, than that 
found at the surface. This lesser modulation would be 
enough to break the symmetry of the BZ and cause a 
visible folding, but not enough to strongly modify FS 
topologies and volumes. The subsurface layer is then 
responsible for all 6-states while the surface for s-states. 
Measured FS volumes [Fig. l(c,d)] support this model: as 
shown in Tab. I, while the s-FS electron counting is some- 
what reduced (n s e ~3.87), the 6-FS electron counting is 
much closer to the dHvA results adding up to ng~3.94, 
i.e. almost 4 electrons as expected in the Ru t<i g or- 
bitals (this still seems to suggest the possibility that these 
near-surface layers are slightly hole-doped, possibly due 
to the cleaving-induced Sr-vacancies revealed by scan- 
ning tunnelling microscopy [20]). Finally photoemission 
from a second layer, also folded but otherwise similar to 
the bulk, could even explain the so-called 'shadow band' 
which remained after high-temperature cleaving had sup- 
pressed the surface states in the original work uniting the 
dHvA and ARPES pictures of the FS of Sr 2 Ru0 4 [5]. 

To conclusively validate this explanation we have per- 
formed a 5-Ru02-layer LDA+SO slab calculation us- 
ing the linearized augmented-plane-wave method in the 
WIEN2K package [21]. Explicitly including these lay- 
ers and the vacuum in the unit cell allows us to study 
both the structural reconstruction as well as effects due 
to broken symmetry in proximity of the surface in a self- 
consistent way. Additionally, this enables us to estimate 
the surface/bulk character of the resulting bands by pro- 
jecting the corresponding wave functions onto the Ru or- 
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Band 


n LDA+SO 


^ DA+bU (eV A) 


ARPES 


Vp (eV A) 




n* HvA 




1.668 ± 0.007 


2.30 ± 0.01 


1.721 ± 0.040 


0.56 ± 0.05 


4.1 ± 0.4 


- 




0.864 ± 0.006 


2.14 ± 0.01 


0.757 ± 0.020 


0.59 ± 0.05 


3.6 ± 0.3 


- 


is, is 


1.340 =b 0.011 


2.18 =b 0.02 


1.396 ± 0.060 


0.42 ± 0.05 


5.2 zb 0.6 


- 


total 


3.872 ± 0.014 




3.874 ± 0.120 








a b ,a b 


1.740 ± 0.003 


2.14 ± 0.02 


1.760 ± 0.040 


1.30 =b 0.60 


1.6 zb 0.8 


1.781 


Pb,& 


0.967 ± 0.022 


2.60 ± 0.25 


0.903 ± 0.020 


0.80 =b 0.10 


3.3 zb 0.5 


0.921 


lb ,7b 


1.252 ± 0.014 


2.21 ± 0.07 


1.280 ± 0.060 


0.76 ± 0.10 


2.9 zb 0.4 


1.346 


total 


3.959 ± 0.026 




3.943 =b 0.120 






4.048 



TABLE I: Carrier concentrations counting electrons, n e — 2 x Afs/Abz with 2 accounting for the spin degeneracy, and Fermi 
velocities along the (0,0) — (tv,tv) direction, vf, as determined by our ARPES and LDA+SO slab calculations for surface and 
bulk-like electronic structure. ARPES-FS volume estimates are from the phenomenological FS in Fig. 1; LDA+SO results were 
obtained for the 24meV shifted chemical potential to match the average electron counting of surface and subsurface FSs as 
determined by ARPES. The dHvA results from Ref. 4, representative of the bulk Luttinger's counting, are also shown. 



bitals for each layer. We used a symmetric slab with 2 
surface layers surrounding 3 bulk-like layers [Fig. 4(a,b)], 
and the surface structure as determined by LEED [14]. 
This consists of a 7.46° rotation and a ~1% elongation 
of the surface Ru02 octahedra, as well as slightly mod- 
ified surface Sr positions. The resulting constant energy 
contour at 24 meV binding energy, corresponding to a FS 
with an averaged electron counting n e — 3.91 to match the 
s-b averaged value experimentally observed by ARPES, 
is shown in Fig. 4(c): 2 sets of red FSs pertaining to the 
2 outer surface layers, and 3 blue FS sets for the 3 bulk 
layers; hybridization between layers does occur at some 
fc-points, where sheets can be seen to change character. 
Although somewhat complicated by hybridization gaps 
appearing at all crossings, and a 5 S pocket of d x 2_ y 2 char- 
acter at (0,0) and (7r,7r) which experimentally is found 




FIG. 4: (color online), (a) Partial structure of the Sr2Ru04 
5-layer slab used for LDA+SO calculations, with correspond- 
ing surface and bulk- like FS. (b) Reconstructed surface of 
Sr2Ru04, with original undistorted unit cell shown solid and 
new shown dashed, (c) FS obtained from the LDA+SO slab 
calculations by shifting the chemical potential 24meV down 
in energy to match the slight averaged hole-doping observed 
by ARPES for Sr 2 Ru04 surface and bulk- like FS. Surface and 
bulk character are shown in red and blue, respectively. 



just above Ep [11], these calculations show good over- 
all agreement with the phenomenological FS concerning 
the splitting, topology, and volume of the various bands 
(Tab. I); as well as - and most importantly - the surface 
versus bulk-like assignment of all detected features. 

To conclude, by means of ARPES and LDA+SO slab 
calculations we have been able to unravel the surface-to- 
bulk progression of the electronic structure in Sr 2 Ru04. 
We find that a (\/2 x<\/2)R45 reconstruction of surface 
and subsurface layers provides a consistent explanation 
of all detected dispersive features in terms of a progres- 
sion of electronic states induced by structural instabili- 
ties, with no evidence for novel phases driven by topo- 
logical bulk properties, ferromagnetic ordering, or the 
interplay between SO and the broken symmetry of the 
surface. This layer-by- layer approach provides the most 
detailed information on FS volumes, Fermi velocities, as 
well as many body renormalizations ^ RPES /^ DA+so , 
for both surface and bulk-like bands (Tab. I). In analogy 
with the recent work on underdoped cuprates [15, 16], 
which found evidence for the enhancement of electronic 
correlations and ordering tendencies in the surface and 
subsurface region, our study of Sr2Ru04 also highlights 
the significantly more correlated character of the top sur- 
face layer bands, which might thus more strongly benefit 
from an LDA+SO+U description [7]. As a matter of fact, 
the inclusion of U = 1 eV (not shown) immediately leads 
to the lifting of the S s FS predicted by LDA+SO, as a 
result of the overall [/-driven bandwidth renormalization. 

We acknowledge C. Bergemann, M.W. Haverkort, R.G. 
Moore, and G.A. Sawatzky for discussions, D. Wong and 
P. Dosanjh for technical assistance. This work was sup- 
ported by the Max Planck - UBC Centre for Quantum 
Materials (A.N.), the CRC and NSERC's Steacie Fellow- 
ship Programs (A.D.), NSERC, CFI, CIFAR Quantum 
Materials, and MEXT KAKENHI (No. 22103002). 



Electronic address: damascelli@physics.ubc.ca 



5 



[1] Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, 

T. Fujita, J. G. Bednorz, and F. Lichtenberg, Nature 

372, 532 (1994). 
[2] Y. Maeno, S. Kittaka, T. Nomura, S. Yonezawa, and 

K. Ishida, Journal of the Physical Society of Japan 81, 

011009 (2012). 

[3] A. P. Mackenzie, S. R. Julian, A. J. Diver, G. J. McMul- 
lan, M. P. Ray, G. G. Lonzarich, Y. Maeno, S. Nishizaki, 
and T. Fujita, Phys. Rev. Lett. 76, 3786 (1996). 

[4] C. Bergemann, S. R. Julian, A. P. Mackenzie, 
S. NishiZaki, and Y. Maeno, Phys. Rev. Lett. 84, 2662 

(2000) . 

[5] A. Damascelli, D. H. Lu, K. M. Shen, N. P. Armitage, 
F. Ronning, D. L. Feng, C. Kim, Z.-X. Shen, T. Kimura, 
Y. Tokura, Z. Q. Mao, and Y. Maeno, Phys. Rev. Lett. 
85, 5194 (2000). 

[6] M. W. Haverkort, I. S. Elfimov, L. H. Tjeng, G. A. 
Sawatzky, and A. Damascelli, Phys. Rev. Lett. 101, 
026406 (2008). 

[7] G.-Q. Liu, V. N. Antonov, O. Jepsen, and O. K. Ander- 
sen., Phys. Rev. Lett. 101, 026408 (2008). 

[8] J. J. Deisz and T. E. Kidd, Phys. Rev. Lett. 107, 277003 
(2011). 

[9] C. M. Puetter and H.-Y. Kee, EPL 98, 27010 (2012). 
[10] S. Kashiwaya, H. Kashiwaya, H. Kambara, T. Furuta, 

H. Yaguchi, Y. Tanaka, and Y. Maeno, Phys. Rev. Lett. 

107, 077003 (2011). 
[11] K. M. Shen, A. Damascelli, D. H. Lu, N. P. Armitage, 

F. Ronning, D. L. Feng, C. Kim, Z.-X. Shen, D. J. 

Singh, I. I. Mazin, S. Nakatsuji, Z. Q. Mao, Y. Maeno, 

T. Kimura, and Y. Tokura, Phys. Rev. B 64, 180502 

(2001) . 

[12] R. Matzdorf, Z. Fang, Ismail, J. Zhang, T. Kimura, 
Y. Tokura, K. Terakura, and E. W. Plummer, Science 
289, 746 (2000). 

[13] V. B. Zabolotnyy, E. Carleschi, T. K. Kim, A. A. Ko- 
rdyuk, J. Trinckauf, J. Geek, D. Evtushinsky, B. P. Doyle, 
R. Fittipaldi, M. Cuoco, A. Vecchione, B. Bchner, and 
S. V. Borisenko, New J. Phys. 14, 063039 (2012). 

[14] R. G. Moore, Manifestations of Broken Symmetry: The 
Surface Phases o/Ca2-xSr x Ru04, Ph.D. thesis, The Uni- 
versity of Tennessee, Knoxville (2006). 

[15] J. A. Rosen, R. Comin, G. Levy, D. Fournier, Z.-H. 
Zhu, B. Ludbrook, C. N. Veenstra, D. Wong, P. Dosanjh, 
Y. Yoshida, H. Eisaki, L. Petaccia, and A. Damascelli, 
ArXiv e-prints (2011), arXiv: 111 1.2673 [cond-mat.supr- 
con] . 

[16] H.-H. Wu, M. Buchholz, C. Trabant, C. Chang, A. Ko- 
marek, F. Heigl, M. Zimmermann, M. Cwik, F. Naka- 
mura, M. Braden, and C. Schler-Langeheine, Nat. Com- 
mun. 3, 1023 (2012). 

[17] T. Ando, T. Akima, Y. Mori, and Y. Maeno, J. Phys. 
Soc. Japan 68, 1651 (1999). 

[18] In analogy with previous work on Bi-cuprates, one 
may wonder if the subsurface folded bands detected in 
Sr2Ru04 might originate from the diffraction of bulk 
photoelectrons traveling through the reconstructed top- 
most surface layer. While the results of surface degrada- 



tion are strong evidence against this, it is nevertheless 
important to explore the parallel with the Bi-cuprates. 
There, two sets of unexpected features are detected in 
addition to the single band of the ideal CuCb square 
plane [22]: these have been referred to in the litera- 
ture as 'shadow bands' and 'umklapp bands' [23]. It is 
the umklapp bands which were initially suspected to be 
due to the diffraction of the CuCb plane photoelectrons 
traveling through the modulated BiO surface layer, since 
they appeared to exhibit multiple positive and negative 
high-order replicas, consistent with a diffraction scenario. 
However, this was later shown not to be the case: all crys- 
tal planes present the same incommensurate modulations 
as the BiO planes, and the umklapp bands therefore stem 
from initial states intrinsic to the modulated Cu02 plane 
and not from extrinsic photoelectron diffraction [15, 24- 
26]. Additionally, it is not with the umklapp but with 
the shadow bands that the subsurface folded bands in 
Sr2Ru04 share their form of reconstruction. Importantly, 
contrary to the umklapp bands, these shadow bands in 
Bi-cuprates were deemed inconsistent with a photoelec- 
tron diffraction scenario; by analogy, the same would hold 
for Sr2Ru04- While initially attributed to antiferromag- 
netic fluctuations giving rise to a (a/2 x v / 2)R45° larger 
unit cell and a correspondingly folded Brillouin zone [22] , 
the shadow bands in Bi-cuprates were conclusively shown 
to arise from the presence of 2 structurally inequivalent 
Cu atoms per orthorhombic unit cell [15, 24-26] . 
[19] L. Walz and F. Lichtenberg, Acta Cryst. C49, 1268 
(1993). 

[20] Y. Pennec, N. Ingle, I. Elfimov, E. Varene, Y. Maeno, 
A. Damascelli, and J. Barth, Phys. Rev. Lett. 101, 
216103 (2008). 

[21] P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, 
and J. Luitz, An Augmented Plane Wave Plus Local Or- 
bitals Program for Calculating Crystal Roperties (Tech- 
nical University of Wien, Vienna, 2001) edited by K. 
Schwarz. 

[22] P. Aebi, J. Osterwalder, P. Schwaller, L. Schlapbach, 
M. Shimoda, T. Mochiku, and K. Kadowaki, Phys. Rev. 
Lett. 72, 2757 (1994). 

[23] A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. 
Phys. 75, 473 (2003). 

[24] A. Mans, I. Santoso, Y. Huang, W. K. Siu, S. Tavaddod, 
V. Arpiainen, M. Lindroos, H. Berger, V. N. Strocov, 
M. Shi, L. Patthey, and M. S. Golden, Phys. Rev. Lett. 
96, 107007 (2006). 

[25] K. Nakayama, T. Sato, T. Dobashi, K. Terashima, 
S. Souma, H. Matsui, T. Takahashi, J. C. Cam- 
puzano, K. Kudo, T. Sasaki, N. Kobayashi, T. Kondo, 
T. Takeuchi, K. Kadowaki, M. Kofu, and K. Hirota, 
Phys. Rev. B 74, 054505 (2006). 

[26] P. D. C. King, J. A. Rosen, W. Meevasana, A. Tamai, 

E. Rozbicki, R. Comin, G. Levy, D. Fournier, Y. Yoshida, 
H. Eisaki, K. M. Shen, N. J. C. Ingle, A. Damascelli, and 

F. Baumberger, Phys. Rev. Lett. 106, 127005 (2011). 



